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ABSTRACT 

We present a near-infrared image of UY Aur, a (y.'9 separated binary system, using 
the Coronagraphic Imager with Adaptive Optics on the Subaru Telescope. Thanks 
to adaptive optics, the spatial resolution of our image was ~ C/.'! in the full width at 
half maximum of the point spread function, the highest achieved. By comparison with 
previous measurements, we estimated that the orbital period is ~ 1640 it 90 yrs and the 
total mass of the binary is ^ 1.73 it 0.29 Mq. The observed i?-band magnitude of the 
secondary varies by as much as 1.3 mag within a decade, while that of the primary is 
rather stable. This inconstancy may arise from photospheric variability caused by an 
uneven accretion rate or from the rotation of the secondary. We detected a half-ring 
shaped circumbinary disk around the binary with a bright southwest part but a barely 
detectable northeast portion. The brightness ratio is > 57 it 5. Its inner radius and 
inclination are about 520 AU and 42° ±3°, respectively. The disk is not uniform but has 
remarkable features, including a clumpy structure along the disk, circumstellar material 
inside the inner cavity, and an extended armlike structure. The circumstellar material 
inside the cavity probably corresponds to a clump or material accreting from the disk 
onto the binary. The armlike structure is a part of the disk, created by the accretion 
from the outer region of the disk or encounters with other stellar systems. 

Subject headings: stars: individual (UY Aur) — stars: pre-main sequence — techniques: 
high angular resolution 
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1. Introduction 

Many T Tauri stars have protoplanetary disks of gas and dust where planets are born. The 
disks have been investigated in various wavelengths. However, the observations have focused mainly 
on the disk around single stars, while more than half of T Tauri stars are binaries (Ghez et al., 
1993; Leinert et al., 1993). Artymowicz & Liibow (1996) predicted that a binary system has two 
kinds of disks: a circumstellar disk associated with each star and a ring-shaped circumbinary disk 
around the binary system. A cavity exists between the circumbinary disk and the central binary. 
A limited number of studies have examined the disks around binary (e.g., GG Tau) systems (Itoh 
et al., 2002; McCabe et al., 2002; Krist et al., 2005). 

UY Aur is a binary system of classical T Tauri stars (CTTSs; Duchene et al., 1999) in the 
Taurus- Auriga star-forming region {d ~ 140 pc; Elias, 1978). The separation between UY Aur A 
(primary) and B (secondary) is 0'.'88 (= 120 AU; Close et al., 1998; hereafter C98). The spectral 
types of UY Aur A and B are estimated to be MO and M2.5, respectively (Hartigan Sz Kenyon, 
2003). 

A circumbinary disk around the UY Aur binary was detected by millimeter, near-infrared, and 
polarimetric observations. Millimeter interferometric observations of ^^CO emission (Duvert et al., 
1998; hereafter D98) revealed that the disk is extended over 20" in diameter; its mass is calculated 
at < 0.01 Mq. D98 assumed that the disk is rotating in Keplerian motion around the binary 
and estimated the total mass of the binary at ~ 1.2 M0. Near-infrared (J-, H-, and K'- bands) 
observations with adaptive optics (C98) resolved the circumbinary disk with an inner cavity. Its 
inclination and inner radius were determined to be ~ 42° and ~ 500 AU, respectively, assuming 
that the inner edge is circular. C98 also detected a clumpy structure along the circumbinary disk, 
which is expected to have a higher dust density than the surrounding area. J-band polarimetric 
observations revealed a centrosymmetric polarization pattern in the circumbinary dust disk (Potter 
et al., 2000). The amplitude of the polarization has a peak at a position angle (PA) = 135°, which is 
consistent with the direction of its semimajor axis detected by C98. Potter et al. (2000) suggested 
that the disk is flat or flared (not spherical) with inhomogeneous dust density. As mentioned above, 
protoplanetary disks around binary systems have peculiar morphology. 

Here we present a high resolution image of the circumbinary disk around the UY Aur binary. 

We describe near-infrared coronagraphic observations in Section 2 and the data reduction proce- 
dures in Section 3. The morphology of the circumbinary disk, as well as astrometry and photometry 
of the binary system, are presented in Section 4. Finally, the results are discussed in Section 5. 

2. Observations 

Near-infrared (i/-band, 1.6 //m) coronagraphic observations of UY Aur were carried out on 
2002 November 23 and 2005 November 9 with the Coronagraphic Imager with Adaptive Optics 
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(CIAO) on the Subaru Telescope (Tamura et al., 2000). A 2^'0 mask in the instrument blocks the 
light of both UY Aur A and B (their apparent separation being ~ 0'.'9). This mask has 2% 
transmission for the central star (Murakawa et al., 2003). A pupil Lyot stop reduces the diffracted 
light. The field of view is about 22" x 22" with a pixel scale of 0'.'02125 pixer^ and 0'.'02133 pixel"^ 
in 2002 and 2005, respectively. 

During the observations in 2005, the sky was clear and the seeing size was ~ 0'.'3. While 
the Pull Width Half Maximum (FWHM) of the point spread function (PSF) of the previous near- 
infrared observations was ~ 0'.'15 (C98), the FWHM of our observations had a spatial resolution of 
0'.'09, the highest achieved, because of adaptive optics (AO). The Strehl ratio was around 0.1. We 
obtained 72 object frames. Integration time was 20 s x 1 coadd for each frame. After obtaining 36 
frames, we performed dithering to remove hot and bad pixels. HD 282615 was observed as a PSF 
reference star before and after the UY Aur observations. We obtained 60 frames of the reference 
star, using a O^'S mask. Integration time was 10 s x 3 coadds for each frame. FS 125 was observed 
as a photometric calibrator (Hawarden et al., 2001). Twilight-flats and dark frames were taken at 
the end of the night. 

During the observations in 2002, the seeing size was ~ 1''0. We obtained 69 object frames, 
with an integration of 20 s x 1 coadd for each frame. We performed dithering after 34 frames were 
obtained. Nine frames of SAO 76648 were obtained as a PSF reference star. Dome-flats were taken 
at the beginning of the night, and dark frames were taken at the end. 

3. Data Reduction 

The object frames were calibrated using the Image Reduction and Analysis Facility (IRAF). 
The frames taken in 2002 were used only for astrometry and photometry because of unstable PSFs 
due to the poor seeing conditions. First, a dark frame was subtracted from the object frames. 
Then the frames were divided by the normalized dome-flat (taken in 2002) or twilight-flat (taken 
in 2005) to remove pixel-to-pixel variation in sensitivity. Hot and bad pixels were removed from 
the divided frames. An average count in the sky region was subtracted. The count was measured 
along a concentric circle (with a radius of ~ 8" and width of ~ l") on the UY Aur binary, except 
for the spider positions. The same reduction procedures were applied for the reference star frames. 

We subtracted the reference star frames from the object frames to detect faint structures 
buried in the halo of the central binary. The peak positions of the binary and the reference star 
were measured with the imexamine task for all frames. The object frames were shifted so that the 
positions of UY Aur A were centered on the image. The reference star frames were duplicated, and 
then each was shifted to adjust the peak positions of the PSFs between the reference star and each 
component of the UY Aur binary. Because the Subaru Telescope is an altazimuth telescope and 
the instrument rotator was used, the PA of the spider changed with time. The shifted reference 
star frames were rotated to adjust the PA of the spider to that of each object frame. The pixel 
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values of the reference star frame were scaled to adjust the flux ratio between each component of 
the UY Aur binary, and the reference star frames were combined. The combined reference star 
frames were renormalized to match the intensities between the reference star frames and the object 
frames at the 0'.'5 x (/.'S region immediately outside the mask northwest of the binary, where the 
effect of the spider is negligible. The reference star frames were subtracted from each object frame. 
Finally, the object frames were combined. It should be noted that the radial profile of the PSF 
outside the mask is independent of both the peak position of the PSF in the mask and of the mask 
size. Artifacts of the mask residuals made by subtracting the PSF reference star are located within 
the mask of the object because the objects and the PSF reference stars are carefully centered in 
the mask. As a consequence, emission structures in the PSF-subtracted image are real structures 
or the residuals of the PSF halo due to imperfections of the AO compensation. 



4. Results 

4.1. Circumbinary disk around the UY Aur binary 

Figure [T] shows the i?-band coronagraphic image of the circumbinary disk around the UY 
Aur binary. This image was not deconvolved using the PSFs of the reference star. A and B 
represent the positions of the primary and the secondary, respectively. The coronagraph and AO 
enabled us to obtain the high spatial resolution (~ O'.'l) image of the disk. If the PSFs of the 
reference star are subtracted too much or too little, the subtracted residual would be a circular 
ring. Because the observed disk has an elliptical shape, it is not an artifact. The disk has five 
remarkably inhomogeneous structures: 

1. A cavity between the binary and the disk. 

2. Significant differences in brightness between the southwest and northeast portions of the disk. 

3. A clumpy structure along the disk ([1] in Fig. [T]). 

4. Circumstellar material in the cavity ([2] in Fig. [T|). 

5. An extended armlike structure ([3] in Fig. [1]). 

C98 fitted an ellipse to the peak intensity of the disk, with a semimajor axis of ~ 500 AU at a 
PA of ~ 135° that may correspond to the inner edge of the disk. If the disk is circular, it is inclined 
by ~ 42°. To derive the disk parameter, we divided the southwest and southeast parts of the disk 
into 10 by 10 pixel regions. Using least-squares fitting of the peak intensity in each region, we drew 
an ellipse with a semimajor axis of ~ 3'.'7 (= 520 AU) at a PA of ^ 135° with an inclination of 
~ 42° ±3°. 
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4.2. Astrometry 

In 2002, the apparent separation and PA between UY Aur A and B were 0'.'886 it O'.'OOl and 
229?13 ± 0?04, respectively. In 2005, they were 0'.'890 ± O'.'OOl and 231?36 ± 0?04, respectively The 
binary system has an inclination of ~ 42° it 3° relative to the line of sight, if it is coplanar with the 
circumbinary disk. As the major axis of the disk is along a PA of ~ 135°, we ascertained that the 
deprojected separation and PA in 2002 were 1'.'191 ± 0'.'02 and 228?08±0?05, respectively In 2005, 
they were 1'.'194 ± 0'.'02 and 229?74 ± 0?05. 

Figure [2] shows the deprojected separation {top) and PA (bottom) of UY Aur A and B. Com- 
parison with previous observations shows that the separation has been constant at 1'.'19 (= 167 
AU) since 1944. UY Aur B has a constant angular velocity of 0?22 it 0?01 yr~^ with respect to A. 
Assuming a circular orbit, the orbital period and total mass, Mx, of the binary are calculated to 
be 1640 ±90 yr (= 360°/0?22 ± 0?01 yr"^) and 1.73 ±0.29 Mq, respectively The calculated Mt is 
nearly consistent with C98's result estimated from the binary motion (1.6llQg7 Mq), but it differs 
from D98's result estimated from the velocity structure of the circumbinary disk (~ 1.2 Mq). The 
assumption of a circular orbit may cause the mass difference. 

4.3. Photometry 

Magnitudes of UY Aur A and B were measured by aperture photometry using the phot task. 
Photometry was performed on the images after the reduction processes, including dark-subtraction, 
flat-fielding by either dome flat or twilight flat, hot and bad pixel rejection, and sky-subtraction. 
We used 0'.'22 and 0'.'18 radius apertures for the 2002 and 2005 images, respectively. The average 
intensity of the region without the circumstellar structures was set as a sky count. The binary 
magnitudes were derived from the measured magnitude under the mask and the Two Micron All- 
Sky Survey (2MASS) magnitude of the PSF reference star. Note that the mask has 2% transmission. 
In 2002, the magnitudes of UY Aur A and B were determined to be 8.11 ± 0.21 mag and 8.54 ± 0.21 
mag, respectively. In 2005, those figures were 7.85 ± 0.10 mag and 9.55 ± 0.13 mag, respectively. 
In 1996, the values were 8.26 ± 0.07 mag and 9.85 ± 0.06 mag, respectively (C98). The -ff-band 
magnitude of UY Aur B varies up to 1.3 mag over a decade, while A is rather stable. 

5. Discussion 

5.1. Asymmetric brightness of the circumbinary disk 

Figure [3] shows the azimuthal intensity profile along the ellipse in Figure [H which is normalized 
by the peak intensity in the circumbinary disk. The southwest side of the disk (225° < PA < 250°) 
is bright, while emission is barely detected toward the northeast portion. The intensity ratio of 
the southwest side to northeast side is > 57 ± 5, while C98 reported that the brightness ratio was 
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2.7. 



We suggest that anisotropic scattering by dust causes the brightness asymmetry of the disk 
and that the southwest side of the disk is brighter due to the forward scattering of the dust. This 
indicates that the southwest side of the disk corresponds to the near side along our line of sight. 

Based on this hypothesis, we described the disk brightness with a modified function of the 
Henyey-Greenstein scattering phase function (McCabe et al., 2002), 



where g is the asymmetry parameter and $(0scat) is the normalized azimuthal intensity of the disk 
along the inner edge as a function of the scattering angle ^scat- The scattering angle ^scat is defined 
as in McCabe et al. (2002), 



where (j)open is the angle related to the disk height, h\ PAq is the position angle of the brightest part 
(= 235?7); i ~ 42° is the disk inclination (see Figured]), j = 1 for cos(PA — PAq) > 0, otherwise 
j = 0. The scattering angle ^scat has minimum at PA = PAq. 

The normalized intensity $(^scat) is presented in Figure O In the case of a geometrically thin 
disk (i.e., (j)open = 0; Fig. (Sj top), the observed intensity profile is not well reproduced with any g. 
In the case of (popen = 45° (i.e., h = 520 AU; Fig. [5l bottom), the observed profile is fitted with 
g 0.9. However, it is hard to imagine a circumbinary disk with such a large asymmetry parameter 
of the particles and such huge disk height. The brightest portion in the circumbinary disk may be 
a local peak of dust distribution. Otherwise, we propose a situation in which the near side is free 
from heavy extinction while the other area is deeply embedded. We cannot exclude the possibility 
that the extinction varies across the disk and induces a variation in the brightness observed across 
the disk. 



The clump is detected along the disk southeast of the binary (PA ~ 156°) and seems to be 
circular. Its diameter is ~ I'.'S (= 180 AU) in the FWHM, which is significantly larger than the 
FWHM of the PSF (~ 14 AU). The apparent distance between the clump and the center of the 
gravity of the binary is ~ 2'.'9 (= 410 AU). We calculated its surface brightness to be 15.3 it 0.1 mag 
arcsec"^, using an aperture diameter of ~ I'.'S and the same sky count as the binary photometry 
(§ 14. 3p . The clump was also resolved by C98, who suggested that it is a local peak of the dust 
density with a calculated diameter of ~ O'.'S and a surface brightness of 16.1 it 0.1 mag arcsec"^ at 
the i?-band. By comparing our image with the ^'^CO-integrated intensity map (Fig. [6]), we found 
that distribution of gas also had a peak near the near-infrared clump. 



$(^scat) oc {l-g^){l + g^ -2gcos9sc^t) 2 



(1) 




5.2. Mass of clump 
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The lower limit of the clump mass was estimated from its luminosity, assuming optically thin 
conditions in near-infrared. Following Lazareff et al. (1990), the clump luminosity Lciump is defined 
as: 

7rr^ f _ _ 

-^ciump = -^^^Lstar / e ^V%e ^""du, (3) 

where Lgtar is the luminosity of the binary, r is the clump radius. D is the deprojected distance 
between the binary and the clump 520 AU), ri is the optical depth between the binary and 
the clump, T2 is the depth between the clump and the observer, p is the lower limit of the clump 
density, is the scattering opacity of gas and dust in the clump, and u is along a light ray emitted 
from the binary. We simplify Eq. (3) with three assumptions. First, the light paths both from the 
binary to the clump and from the clump to the observer are optically thin (i.e. e~'^^ and e~'^^ are ~ 
1). Second, the clump is spherical, with r ~ 90 AU. Assuming uniform distribution of the dust in 
the clump and single scattering, the average path length of light in the clump is nearly equal to its 
diameter, 2r ~ 180 AU. Finally, because the distance between the binary and the observer is quite 
large (~ 140 pc), luminosity and flux are related as 

-^clump -^clump 

With these assumptions Eq. (3) is transformed into 

-^clump 2Z? , , 

We derived -Fdump = 5.9 x 10^ counts and Kjtar = 4.9 x 10^ counts from the photometry above. 
Given that Ks = 18.5 cm^ g~^ for the mixture of graphite + silicate grains and gas (Wolfire &; 
Cassinelli, 1986), p is calculated to be ~ 5.4 x 10~^® g cm~^ and the lower limit of the clump mass 
is ~ 9.6 M®. High spatial resolution (sub-)mm wavelength observations are required to obtain an 
accurate mass for the clump. 



5.3. Circumstellar material 

Possible emissions were detected inside the cavity in the northwest and southeast areas of the 
UY Aur binary (hereafter called NW and SE, respectively). They may be spiral arms extending 
from the circumbinary disk to the central binary. C98 reported the presence of an arm at O'.'S 
from UY Aur B (see the J-band image in their Fig. 1). NW nearly corresponds to the position 
of C98's arm. Artymowicz & Lubow (1996) demonstrated that the secondary star perturbing a 
circumbinary disk produces such spiral arms. Otherwise, SE could be a clump. We calculated its 
lower mass limit to be ~ 6 x 10~^ M^. Data reductions produced negative artifacts as strong as 
positive ones in the region just outside the mask (—30° < PA < 90°), so discussions about NW and 
SE are less reliable. 
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5.4. Armlike structure 

The armlike structure was detected at PA ~ 270° (see [3] in Figure [1]), which is elongated to 
the north and extends to ~ 6'.'5 (= 910 AU) from the UY Aur binary. Given that the structure is 
a part of the circumbinary disk, we considered two possible mechanisms for its origin. First, cold 
dust may have accreted from the outer region of the circumbinary disk to the central binary. The 
outer radius of the disk was ascertained to be > 10" (= 1400 AU, D98), with the circumbinary 
disk having higher gas surface density than the other area (see Fig. 2 in Ochi et al., 2005). Second, 
the structure may have formed by an encounter with another stellar system. Pfalzner (2003) 
demonstrated that a close flyby of a star creates an arm structure in a circumstellar disk. We 
investigated the Taurus-Auriga star-forming region. Our calculation, which considered the proper 
motion of UY Aur and other young stellar objects (Jones & Herbig, 1979), showed that the paths 
of RW Aur and JH 218 intersected with that of UY Aur during the past 10^ yrs. However, none of 
these stars show evidence of a stellar encounter, such as an armed disk structure. 

5.5. Light variation of UY Aur B 

Brandeker et al. (2003) reported that the -fT-band flux ratio of UY Aur A to B changes from 
4.33 ± 0.36 (1996 Oct; C98) to 1.47 ± 0.01 (2000 Feb; Brandeker et al., 2003). UY Aur B had 
a magnitude variation up to 1.3 mag at the i?-band during the decade, while the magnitude of 
A was rather stable. This variation probably originates from the photospheric variability of UY 
Aur B. The light variation often observed in CTTSs may be caused by variable accretion rates to 
the central star or by rotation of a star with a hot spot on its surface. Simultaneous multiband 
observations will clarify the source of this variation. 

6. Conclusions 

Using the Subaru/CIAO, we obtained the highest spatial resolution (~ O'.'l) near-infrared coro- 
nagraphic observations yet achieved of the UY Aur binary. The southwest side of the circumbinary 
disk around the binary was clearly seen, while the northeast side was barely detected. Its inner 
radius and inclination are ~ 520 AU and ~ 42°, respectively. The following disk features were 
observed: 

1. The intensity ratio (near/far) of the inner disk edge is calculated to be > 57 ± 5. We de- 
termined that the southwest bright side of the disk is located near to us, assuming forward 
scattering of the dust. By fitting the Henyey-Greenstein phase function given a disk height 
of ~ 520 AU, we derived the asymmetry parameter g ~ 0.9. Inhomogeneous distributions of 
dust in the disk and in front of the disk may also account for the huge ratio. 

2. A clumpy structure was resolved in the southeast region of the disk. The apparent distance 
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from the binary was ~ 410 AU. This is probably a dense clump of dust and gas. The lower 
limit of the clump mass is ~ 9.6 M®, assuming optically thin conditions. 

3. Circumstellar materials were detected at two regions within the inner cavity. The material in 
one region may have accreted from the circumbinary disk to the central binary through the 
cavity. The second region may also be a clump. However, these structures are located just 
outside the coronagraphic mask and further observations are needed. 

4. An armlike structure northwest of the binary was probably created by dust accretion from 
the outer region of the disk or through a stellar encounter. We detected two candidates with 
proper motions that would have allowed them to have encountered the binary. 

We also conducted astrometry and photometry of the binary. The deprojected separation 
between UY Aur A and B was constant at ~ 1'.'19 (= 167 AU) over about 60 yrs, suggesting 
that the binary has a circular orbit. The total mass and the orbital period of the binary were 
~ 1.73 ±0.29 Mq and ~ 1640 yrs, respectively. The iJ-band magnitude variation of UY Aur B was 
up to 1.3 mag over a decade. This variation may be caused by UY Aur B's photospheric variability. 
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Fig. 1. — i/-band coronagraphic image of the circumbinary disk around the UY Aur binary. The 
field of view is 13'.'2 x 13'.'2. Diameter of the occulting mask is 2'.'0 (= 280 AU). The PSFs of the 
central binary are subtracted. The figure on the right side is the original image (left) with marks. 
A and B in the mask represent the positions of the UY Aur primary and secondary, respectively. 
The southwest side of the disk is clearly seen while the northeast side is barely detected. The 
southwest side corresponds to the side nearest to us. The solid ellipse indicates the inner edge 
of the disk (r ~ 520 AU). The disk inclination i was calculated to be ~ 42° it 3°. [1]: a clump, 
the brightest portion along the circumbinary disk; [2]: circumstellar material in the cavity; [3]: an 
extended armlike structure. 
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Fig. 2. — Top: Deprojected separation between UY Aur A and B. If a circular orbit is assumed, its 
radius is ~ 1'.'19 (= 167 AU) (described by the dotted line). Data are from Joy & Van Biesbroeck, 
1944; Ghez et al., 1993; Leinert et al., 1993; Ghez et al., 1995; C98; White & Ghez, 2001; Hartigan 
&; Kenyon, 2003; McCabe et al., 2006; Brandeker et al., 2003; and this work (last two points). 
Bottom: Deprojected position angle of the binary. Data are from the same sources as in Figure [2] 
(top). The dotted line drawn by least-square fit indicates that UY Aur B has an angular velocity 
of ~ 0?22 yr^^ with respect to A. 
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Fig. 3. — Azimuthal intensity profile of tlie inner edge of the circumbinary disk. The points 
represent the average intensity of the elhpse (Fig. 1) in 3 x 3 pixel aperture. The intensity 
is normalized by the peak intensity (PA ~ 156°) in the disk, which corresponds to the clump. 
Forward scattering of the dust creates the peak in the near side (225° < PA < 250°). In contrast, 
the far side (45° < PA < 70°) is barely detected. The intensity ratio (near/far) is calculated to be 
> 57 ± 5. The vertical dotted lines represent the position angles of the spider. North is PA = 0°, 
and east is PA = 90°. 
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Fig. 4. — Schematic view of the model of the UY Aur binary and its circumbinary disk. The 
scattering angle ^scat, disk height h, and <?!)open are shown. 
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Fig. 5. — The observed azimuthal intensity profile of the circumbinary disk {points) as a function 
of scattering angle ^scat- Lines represent the Henyey-Greenstein function (dashed lines for = 0.1, 
dotted lines for g = 0.5, and solid lines for g = 0.9). The small scattering angle corresponds to 

the near side of the disk while the large angle corresponds to the far side. The local peaks at ^scat 
= 65°-95° (top) and 20°-50° (bottom) represent the clump. Top: Flat disk (h = 0). The function 
does not fit well to the observed profile with any g. Bottom: Geometrically thick disk {h = 520 
AU). The observed profile is fitted by the Henyey-Greenstein function with g = 0.9. 



Fig. 6. — //-band image overlaid ^^CO (J = 2 ^ 1) integrated intensity map (D98; contour). The 
spatial resolution in ^^CO (J = 2 ^ 1) is 2'.'9 x 2'.'3. The field of view is 14'.'4 x 14'.'4. The gas 
distribution has a peak near the near-infrared clump. ^^CO emission extends over 20" (= 2800 
AU) in diameter. 



